A new treatment of the in-medium chiral condensates 
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A new formalism to calculate the in-medium chiral condensate is presented. At lower densities, this approach 
leads to a linear expression. If we demand a compatibility with the famous model-independent result, then the 
pion-nucleon sigma term should be six times the average current mass of light quarks. QCD-like interactions may 
slow the decreasing behavior of the condensate with increasing densities, compared with the linear extrapolation, 
if densities are lower than twice the nuclear saturation density. At higher densities, the condensate vanishes 
inevitably. 



The behavior of chiral condensates in a medium 
has been an interesting topic in nuclear physics 
[1] . A popular method to calculate the in- medium 
quark condensate is the Feynman-Helmann the- 
orem. The main difficulty, however, is the as- 
sumptions we have to make on the derivatives of 
model parameters with respect to the quark cur- 
rent mass. 

To bypass this difficulty, we will apply a similar 
idea as in the study of strange quark matter [2- 
5] by defining an equivalent mass. A differential 
equation which determines the equivalent mass 
will be derived. At lower densities, the new for- 
malism leads to a linear decreasing condensate. 
A comparison with the result in nuclear matter 
implies that the pion-nucleon sigma term should 
be six times the average current mass of light 
quarks. At higher densities, it turns out that the 
decreasing speed of the condensate with increas- 
ing densities is lowered, compared with the linear 
extrapolation. 

The QCD Hamiltonian density can be schemat- 
ically written as 

Hqcb ^ H]^ + '^'rnioqiqi + Hi, (1) 

i 

where Hy^ is the kinetic term, m^o is the current 
mass of quark flavor i, and Hi is the interacting 
part of the Hamiltonian. The sum goes over all 
flavors involved. 



The basic idea of the mass-density-dependent 
model of quark matter is that the system energy 
can be expressed as in a noninteracting system, 
where the strong interaction implies a variation 
of the quark masses with density. In order not to 
confuse with other mass concepts, let us call such 
a density-dependent mass as an equivalent mass. 
It can be separated into two parts, i.e., 

= rriiQ -f mj, (2) 

where the first term is the quark current mass and 
the second part is a flavor independent interact- 
ing part. Therefore, we will have a Hamiltonian 
density of the form 

-ffcqv = -ffk + ^ miqiqi. (3) 

i 

We require that the two Hamiltonian densities 
i?oqv and Hqcd should have the same expecta- 
tion value for any state j'l'), i.e., 

(^liJcqvl^-) = {^\Hcicd\^). (4) 

Applying this equality to the state jne) with 
baryon number density ub and to the vacuum 
state |0), we have 

{Hcqv)nB ^ {Hcqv)o = (^^QCD)riB ^ (-H^QCd)o- (5) 

Here we use {A)n-g={nB\A\nB) and (^)o=(0|^|0) 
for an arbitrary operator A. 
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We restrict ourselves to systems with uni- 
formly distributed particles where we can write 
{^\m{nB)qq\'^) = m{nB){'^\qq\'^)- Accordingly 
we can solve Eq. (5) for mi, getting 

"^i = ^//- \ — /- V X . (6) 

22i{{qiqi)nB - {qiqi)o) 

where ej = {Hi)n^ — (-ffi)o is the interacting en- 
ergy density. 

Therefore, considering the quarks as a free sys- 
tem, while keeping the system energy unchanged, 
they should acquire an equivalent mass corre- 
sponding to the current mass plus the common 
interacting part shown in Eq. (6). The equiva- 
lent mass is a function of the quark current mass 
and of the density. Note that quark confinement 
implies the following natural requirement: 

lim mi = 00. (7) 

Because the Hamiltonian density ileqv has the 

form of a system of free particles with equivalent 
masses mi, the energy density of quark matter 
can be expressed as 

e = > I \/ k'^ + m'? k'^dk 



(8) 



where g = 3 (colors) X 2 (spins) = 6 is the degener- 
acy factor, and 

is the Fermi momentum of the quarks. The func- 
tion F{x) is 



F{x) = %Vx^ + 1 (2a;2 + 1) - sh.-\x)]/x^.{W) 
8 

For convenience, let us define the function f{x) 

f{x) = -x'^d[F{x)/x]/dx 

= ^ [.Tv/a;2 + l-sh-i(x)] /x\ (11) 

On the other hand, the total energy density can 
also be expressed as 

e = y ^1 \/k'^ + m%k'^dk + ei 



3 \- TP f 

— riB > ruioF 

N{ ^ \mio 



(12) 



where iVf = 2 is the number of flavors, the first 
term is the energy density without interactions, 
and the second term is the interacting part. Be- 
cause ns/Nf is the baryon number density for 
each flavor, the Fermi momentum fcfo here is 



, 18 2"-B 
kiO = — TT — 

g Ni 



1/3 



(13) 



It looks similar to the non-interacting system. 
However, the Fermi momentum fcf in Eq. (9) is 
bigger. It has been boosted because of the Fermi 
momentum dependence on density through the 
equivalent mass. In the appendix A, we will give 
a proof for the boosting of the Fermi momentum. 
Combining Eqs. (12) and (8) we identify 



ei 
3nB 



E 



niiF 



The Hellmann-Feynman theorem gives 
(f|^i7(A)|vI/) = |-(*|i?(A)|vI/), 



(14) 



(15) 



where j^f') is a normalized eigenvector of the 
Hamiltonian H[X) which depends on a parameter 
A. 

On application, in Eq. (15), of the substitu- 
tions A TOio and H{X) J d^xHQCB, one 
gets (*|/d3a;g-i(z,|*) = / d^a^FqcDl*) 

for each flavor i. Applying this equality, re- 
spectively, to the state Iub) (quark matter with 
baryon number density ne) and to the vacuum 
|0), one obtains 



{qiqi)nB - {qiqi)o = 



de 
dmio ' 



(16) 



where e = (i?QCD)nB - {Hqgd)o is the total en- 
ergy density. Now let us substitute Eq. (8) into 
Eq. (16), carry out the corresponding derivative, 
and sum over the flavor index. We get 



i 



(17) 
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with V = X]j5/9mjo. Note that V is a diflFeren- 
tial operator in mass space. 

Comparing this equation with Eq. (6) we have 



Vmi 



ei/(3nB) 



1. 



(18) 



Replacing ei/(3riB here by the right hand side 
of Eq. (14) we get a first order differential equa- 
tion for the interacting equivalent mass. 

Such a mass really exists, and we can prove that 
it can be expressed in terms of the interacting 
energy density ei formally as 



mi 



ei/(3nB) 



3nB ' A^i 



Vf Ei / 



(19) 



In the flavor symmetric case, i.e., muo = rudo = 
■■■ = mo, we have ruu = rud = ■ ■ ■ = m, {quQu) = 
{OdQd) = ■ ■ ■ = {qq), and V = d/dmo. In this case 



dm__rnF (fcf/m) - {ho/mo) 
dmo 



mif{k{/m) 
1 ei 



1, 



{qq)o N{{qq)o mi 



(20) 



(21) 



(22) 



Since limj-^o -^(a;) = 1, Eq. (22) becomes at 
lower densities 



m = mo + 



ei 



SNfUB 



(23) 



This means mi = ei/(3A'^fnB), i.e., ei/mi = 
3A'fnB. Substituting this ratio into Eq. (21), we 
get 



{qq)r 



{qq)o 

with 

* 

n = - 



^ riB 



6mo 



-3(99)0 



(24) 



(25) 



where m^^ k, 140 MeV is the pion mass and w 
93.2 MeV is the pion decay constant. 

Since we have said nothing about the form of 
the interacting energy density, our result is model 



independent. Recalling that there is a model- 
dependent result in nuclear matter, i.e., 



Mii = i-^withp* 

{q<i)a P* 



(26) 



first proposed by Drukarev et al. [6], and later 
re-justified by many authors [7], we get, from the 
requirement n* = p* , the very interesting relation 
itn = 6mo, i.e., the pion-nucleon sigma term un 
is six times the average current quark mass mo. If 
one takes ctn = 45 MeV [8-10] and mo = (m„o + 
TOdo)/2 = (5 + 10)/2 = 7.5 MeV [11], we confirm 
this result. 

The chiral condensate at higher densities can 
be calculated from Eqs. (20)-(22) if we know 
the interacting energy density ej from a realis- 
tic quark model. In the following, we consider a 
simple example. 

Denoting the average distance between quarks 
by f, the interaction between quarks by 
v(mo,nB), and assuming that each quark can 
only interact strongly with other A^o nearest 
quarks at any moment, because of the satiiration 
of strong interactions, the interacting energy den- 
sity ej can be linked to density by 



ei = -A^o'^Bv(mo,r). 



(27) 



The average inter-quark distance r is linked to 

1 /3 

density through f = ^/rif, . Here ^ is a geomet- 
rical factor related to the way in which we group 
the quarks together. In what follows, we have di- 
vided the system into sub cubic boxes, being then 
^ = 1/3^/^. We will take Nq = 2 since a quark 
has a trend to interact strongly with other two 
quarks to form a baryon. The concrete value of 
A^o as well as the value of ^ have only a marginal 
influence on the density behavior of the chiral 
condensate. 

Substituting Eq. (27) into Eqs. (21) and (22), 
we have, respectively. 



(gg)«B 
{qq)o 



1 



Nq Ub v 

2N{ n* mi ' 



mF 





-^f( 











No 
2Nf 



(28) 



v(mo,nB). (29) 
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If the parameter Nq diverges faster than fcf or 



1/3 



at extremely higher densities, we have 



hm v(mo, r) = 0. 



(30) 



which is consistent with asymptotic freedom. 

To solve Eq. (20), we need an initial condition 
at mo = mg. Let us suppose it to be 



m{mo,nB) = m(nB). 
Usually, we will have 
v(mo,nB)Uo=m* = v(f). 



(31) 



(32) 



where v(f) is the inter-quark interaction for the 
special value ttiq of the quark current mass niQ. 

Eq. (20) is difficult to solve analytically. How- 
ever, this can be done at lower densities. 

Let's rewrite Eq. (19) as 



mi 



v(too, 



_2_ f /fetfl.'i -L gV("'o.'') 
No-' \mo J 



(33) 



dmo 



At lower densities, the Fermi momentum k{ 
is small, so the function F{x) approaches to 
1. Accordingly, from Eq. (29) we get mi = 
^^v(mo, tt-b). Replacing the left hand side of Eq. 
(33) with this expression, and integrating the re- 
sulting equation under the initial condition given 
in Eq. (32), we have 



mi(mo,nB) 



1 



Nf'' I mo 



dmo. (34) 



In general, an explicit analytical solution for 
the condensate is not available, and we have to 
perform numerical calculations. For a given inter- 
quark interaction v(f), we can first solve Eq. (29) 
to obtain the initial condition in Eq. (31) for the 
equivalent mass, then solve the differential Eq. 
(20), and finally calculate the quark condensate 
through Eq. (28). 

There are various expressions for v(f ) in litera- 
ture, e.g., the Cornell potential [12], the Richard- 
son potential [13], the so-called QCD potentials 
[14,15], etc. They are all flavor- independent. 
Let's take a QCD-like interaction of the form 



v(r) 



_ 4as(r) 

ar 

3 f 



(35) 



The first term af is the long-range confining part. 
The second term incorporates perturbative ef- 
fects. To second order in perturbation theory, 
one has [14,15] 



boX{r) 



A(f) A(f) 



(36) 



where [16] 

A(f) = ln[(fA^)-2 + 6] 



(37) 



and6o = (llAc-2Aff)/3, h = [34:N^-N{{13N^- 
3)/Afc]/3, and 62 = (31iVc - 10iVf)/(96o) +27E for 
SU(A/c) and A^f flavors. 7e is the Euler constant. 

Besides these constants, there are three free pa- 
rameters, i.e. a, Atos, and b. The QCD scale pa- 
rameter is usually taken to be Aj^j = 300 MeV. 
The value for the string tension a from potential 
models varies in the range 0.18 — 0.22 GeV^ [17], 
and we take a = 0.2 GeV^. As for the param- 
eter b, we take three values i.e. 10, 20, and 30, 
in the reasonable range [16]. The value of ttiq in 
Eq. (32) is taken to be 7.5 MeV. The numerical 
results are plotted in Fig. 1. 

In Fig. 1, the straight line is the linear extrap- 
olation of Eq. (24). It does not depend on the 
form of the inter-quark interaction v(f), and so 
is 'model-independent'. The other three lines are 
for Too = 7.5 McV, but for different b values. At 
lower densities, the chiral condensate decreases 
linearly with increasing densities. When the den- 
sity increases, being less than two times the nu- 
clear saturation density, the decreasing speed is 
slowed. However, for even higher densities, it can 
be shown that the condensate vanishes rapidly. 

It should be noted that if the Fermi momen- 
tum in Eqs. (8) had not been boosted, the main 
Eq. (21) is stm valid while Eqs. (20) and (22), 
and accordingly Eqs. (23), (29), (33)~(34), and 
especially the important equation (24) and the re- 
lation between the pion-nuclcon signia term and 
the current quark masses would be different by a 
factor, as has been formulated in the first part of 
Ref. [18]. 
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Figure 1. Density dependence of the quark con- 
densate in quark matter. 



A. Why should the effective Fermi mo- 
mentum be boosted 

In this appendix, we show that the effective 

Fermi momentum in the equivalent mass ap- 
proach should be boosted to a higher value. 
We start from 



d{VE) = Td{VS) - PdV + ^id{Vn), 



(38) 



which is the combination of the first and second 
laws of thermodynamics. Here n is the particle 
number density, E is the energy density, and S is 
the entropy density. Because the system is uni- 
formly distributed, the corresponding extensive 
quantities are, respectively, Vn, VE, and VS. ji 
is the chemical potential. From this expression 
we can get 



T = 



dE 



dS „ 
P + E-TS-nn 



-V 



dE 
dV 



(39) 
0, (40) 



dE 



(41) 



At zero temperature, the entropy becomes zero, 
Eqs. (40) and (41) become, respectively, 

P = -E + nn, (42) 
dn = dE/ii. (43) 

In our equivalent mass model, the energy den- 
sity is given with the equivalent mass as 

E = £^ r + m-'p^dx 



(44) 



167r2 



kfJkf 



to2 {2kf 



m 



-mSh-i 



(45) 



where the Fermi momentum kf satisfies fcf 

\JliF--rn?- . The degeneracy factor g* 



IS 



A^f (flavor) X 2 (spin) X 3 (color) . 

iProm Eq. (45) we get dE = ||dfcf + ||dm. 
Substituting this into Eq. (43) then gives 

dn 

rn^ sh^ ^ (ki/ni) 



+ 



g*m 



47r2 



kf 



dm. (46) 



If the mass does not depend on the density or 
Fermi momentum, the second term vanishes. One 
then has 



ki = 



9* 



1/3 



(47) 



Eq. (47) is the well-known expression for the 
non-interacting system. However, in the mass- 
density-dcpendcnt case where interactions are 
treated non-pertutbatively by defining an equiv- 
alent mass, the quark number density should be 
given by integrating over both sides of Eq. (46): 



n 



■ 



9*kf 
67r2 

47r2y 



kf 



m^sh ^(fcf/m) 



mdm (48) 



S,n 



Usually the equivalent mass is a big quantity, 
much larger than the current mass. Therefore, 
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the ratio fcf/m is small if the densities are not too 
high. Let us then expand the integrand of the 
second term on the right hand side of Eq. (48) 
with respect to kf/m, taking then the lowest or- 
der term. We get: 



9* ,3 9* f , 



(49) 



Because of quark confinement and asymptotic 
freedom, m increases with decreasing k[. There- 
fore, the simplest parametrization should be 

C 

TO = Too + — g-. (50) 

with C being a constant. To be consistent with 
the linear confinement, the exponent Z is equal to 
1 ^. However, to reproduce the presently accepted 
value for the pion-nucleon term (about 45 MeV), 
Z should be about 3/2. Substituting Eq. (50) into 
Eq. (49) then gives 

1/3 



(51) 



Comparing Eqs. (47) and (51), it is obvious 
that, for the same density, the Fermi momentum 
of the interacting system is different from that of 
the non-interacting case. When taking Z = 3/2, 
Eq. (51) becomes Eq. (9). 
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